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Abstract 

It has been known for some time that the well established LMA solution to 

the observed solar neutrino deficit fails to predict a flat energy spectrum for Su- 

PLh' perKamiokande as opposed to what the data indicates. It also leads to a Chlorine rate 
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which appears to be too high as compared to the data. We investigate the possible 



D I solution to these inconsistencies with non standard neutrino interactions, assuming 



that they come as extra contributions to the f^i^^ and UaC vertices that affect both 
the propagation of neutrinos in the sun and their detection. We find that, among 
the many possibilities for non standard couplings, only the diagonal imaginary ones 
t~-. ' lead to a solution to the tension between the LMA predictions and the data, implying 

^«0 . neutrino instability in the solar matter. Unitarity requirements further restrict the 

solution and a neutrino decay into an antineutrino and a majoron within the sun is 
the one favoured. Antineutrino probability is however too small to open the possibility 



f~^ ■ of experimentally observing antineutrinos from the sun due to NSI. 
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1 Introduction 

Neutrino non-standard interactions (NSI) have been introduced long ago [Tl[2] to account 
for a possible alternative solution to the solar neutrino problem. Since then a great deal 
of effort has been dedicated to study its possible consequences. To this end possible NSI 
signatures in neutrino processes have been investigated, models for neutrino NSI have been 
developed and bounds have been derived [3HT7]. Specific investigations of NSI in matter 
have also been performed within the context of supernova ^^ and solar neutrinos [T9H23] . 

Although LMA is generally accepted as the dominant solution to the solar neutrino 
problem ^25i;26j, not only its robustness has been challenged by NSI, as it can shift the 
LMA solution to the dark side region of parameter space [22] , but also some inconsistencies 
remain regarding its agreement with the data [271128] . In fact, while the SuperKamiokande 
(SK) energy spectrum appears to be flat [29|[30] . the LMA prediction shows a clear negative 
slope in the same energy range. With the expected improvement in the trigger efficiency 
for threshold electron energies as low as 3 MeV to be reached in the near future [31], such 
a disagreement, if it persists, may become critical. Moreover the LMA solution predicts an 
event rate for the CI experiment [32] which is 2a above the observed one [32] • These are 
motivations to consider 'beyond LMA' solar neutrino solutions in which NSI may play a 
subdominant, although important role. 

In order for NSI to be detectable and therefore relevant in physical processes, the char- 
acteristic scale of the new physics must not be too much higher than the scale of the physics 
giving rise to the Standard Model interactions, Aew — Gp ■ Possible realisations are 
one loop radiative models of Majorana neutrino mass [31], supersymmetric SO(IO) with 
broken D-parity [35], the inverse seesaw in a supersymmetry context [36] or triplet see- 
saw models [37]. Since the scale at which the new interaction arises is supposed to be 
not too far from the electroweak scale, its coupling may be parameterised by Gpe where 
e ~ Al^/A^p ~ 10^^ for d=6 or £ ~ A|;^/Aj^p ~ 10""^ for d=8 operators respectively. 
For type I seesaw, NSI are of course negligible. 

In this paper we will be concerned with NSI both at the level of propagation through 
solar matter and at the level of detection. Matter NSI are defined through the addition of 
an effective operator to the Lagrangian density 

C'^^si = -^^Gpei^^ifYPf] [vcl.PLVp] (1) 

where / = e,u,d, P denotes the projection operator for left and right chirality and ej^i^ 
parameterizes the deviation from the standard interactions. At present there is no evidence 
at all of such operators generated at a scale A^p, hence the variety of theoretical models 
for the physics accessible to the LHC 

There are several ways to introduce NSI. For instance in fermionic seesaw models, once 
the heavy fermions (singlets or triplets) are integrated out, modified couplings of leptons to 
gauge bosons are obtained in the form of a non-unitary leptonic matrix. The strong bounds 
on the deviation from the unitarity of this matrix constrain these NSI to be < O(10~^) 
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Alternatively NSI can be generated by other new physics above the electroweak scale not 
related to neutrino masses. As a consequence an SU(2) gauge invariant formulation of NSI 
is required, since any gauge theory beyond the standard model must necessarily respect its 
gauge symmetry. Strong bounds from four charged fermionic processes [38] and electroweak 
precision tests requiring fine tunings imply that possibilities are limited for such scenarios 
[111139]. Another way to introduce NSI is by assuming that there are extra contributions to 
the vertices z/q,i/^ and z/^e. In such a case the parameters e^/j describe the deviation from the 
standard model vertices and are treated like the standard interactions. It is possible that 
other effects are present in this case that depend on the nature and number of particles that 
may be introduced in a particular model. We adopt this procedure in the present paper and 
assume these model dependent effects to be negligible. 

The paper is organized as follows: section [2] is devoted to the study of the propagation 
and detection of solar neutrinos. We start by reviewing the derivation of the neutrino 
refraction indexes with standard interactions (SI) and their generalization to NSI in order 
to obtain the matter Hamiltonian. The survival and conversion probabilities to z/^ and 
v-r are then evaluated through the numerical integration of the Schrodinger like equation 
using the Runge-Kutta method and the experimental event rates are obtained. We use the 
reference solar model with high metalicity, BPS08(GS) [ID]. In section [3] we investigate 
the influence of the NSI couplings on these rates in order to find whether and how the 
fits can be improved with respect to the LMA ones. We concentrate in particular on the 
elimination of the upturn in the SuperKamiokande spectrum predicted by LMA for energies 
below 8-10 MeV not supported by the data [271430] and on the Chlorine rate whose LMA 
prediction exceeds the data by 2a [32l|33] . We find the intriguing result that only imaginary 
diagonal couplings in the matter Hamiltonian provide a suitable solution to the problem. 
In fact all other possibilities, namely imaginary off diagonal couplings or real couplings, 
either diagonal or off diagonal, are unable to change the LMA solution. This means that 
eliminating the tension with the data implies complex matter Hamiltonian eigenvalues and 
therefore neutrino decay induced by the solar matter. In section H] it will be seen that 
the condition of consistency with the data on the event rates is not sufficient, owing to 
the requirement of unitarity which severely restricts the correct solution. Furthermore, the 
neutrino decay in the sun suggested long ago into an antineutrino and a scalar particle 
(majoron) [HHS] is the favoured channel. In section [S] we present a brief discussion and 
summarize our main conclusions. 



2 Interaction potentials, the Hamiltonian and the rates 

In this section we develop the framework that will be used as the starting point for the anal- 
ysis of the NSI couplings in section [31 To this end we review the derivation of the neutrino 
interaction potentials in solar matter, its generalization to non standard interactions along 
with the corresponding matter Hamiltonian and the event rates. 



2.1 Interaction potentials and the Hamiltonian 

While UgS propagate through solar matter their interaction with electrons proceeds both 
through charged and neutral currents (CC) and (NC). Recalling that for standard interac- 
tions (SI) each tree level vertex accounts for a factor 

^^ iT,L-2Qjstn^ew), (2) 
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inserting the W, Z propagators and the electron external lines, one gets for the fg interaction 
potentials 

{Ve)cc = GfV2N, , {Ve)NC = G p j ^{-\ + ^sir?e^)N, (3) 

where Gp is the Fermi constant, Gpj ^2 = g'j^/Sm^. For the interactions with quarks 
only neutral currents are involved and the additivity of the quark-current vertices gives for 
protons 

Vp = {V^)nc = Gf/V2{1 - Asin^ew)N, (4) 

and for neutrons 

Vn = {yn)NC = -GF/V2Nn. (5) 

Hence the neutrino interaction potential is for standard interaction^ 



ViSI) = V, + Vp + Vn = GfV2N, (l - ^\ =V, + Vn 



(6) 



with K = (K)cc + (K)7vc and V, = V, + V^ = GpV^N,. 

In order to introduce NSI we assume that each diagram associated to neutrino propaga- 
tion in matter (i.e. CC and NC currents in Va e~ and NC currents in z/^ m, z/^ d scattering) 
is multiplied by a factor e^'^' parameterising the deviation from the standard model. So 
we assume that the interaction potential for v^ (« = e, /x, r) on electrons involves both CC 
and NC giving rise to possible lepton flavour violation: z/^ for a ^ e may have CC. So for 
the charged current of p^, with electrons we have 

{Ve)cc{NSI) = ^{e%)ccK (7) 
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and for the neutral current 



{ye)Nc{NSI) = /V(-l + ^sin^Ow){e%)NcNe (8) 

where the NSI couplings affecting the CC and NC processes should in principle be distin- 
guished. 



^AU expressions are divided by 2 to account for the fact that the medium is unpolarized. 



Using equation (|2]) and the additivity of the quark-current vertices, one gets for the 
neutrino interaction potential with protons 



Vp{NSI) 
Similarly for neutrons 
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In both (|9]) and ( ITOl) only neutral currents are involved. 

Adding (ED, (E]), (E]) and ((lO]) and dividing by 2 one finally gets 
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V{NSI) = GFV2Ne 
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SO that the full interaction potential is the sum of eqs.(l6]) and ( ITT1) . 

In the case of the standard interactions, the interaction potentials for z/e and Va constitute 
a diagonal matrix because they cannot be responsible for flavour change [eq.([6])]. This may 
occur as a consequence of the vacuum mixing angle (oscillations) [25, 26] or the magnetic 
moment for instance p7lf28] . On the other hand, in the case of NSI the interaction potentials 
[eq.l fTT]) ] constitute a full matrix in neutrino flavour space. 

In order to obtain the matter Hamiltonian eqs.(l6]) and (ITTi) must now be added. In the 
flavour basis this is 

1 \ / Vee{NSI) 

'Hm = V,\ I + I V,, 
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where in the first term, describing the standard interactions, the additive quantity Vn which 

is proportional to the identity, has been removed from the diagonal. In the second term Vap 

{a, /3 = e, /i, r) denote the matrix elements of the interaction potential matrix flTTl) . Finally 

in the mass basis 
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where U is the PMNS matrix 

nii {i = 1,2,3) the neutrino mass. Upon insertion of this Hamiltonian expression in the 



m"^ with 



neutrino evolution equation, the survival (-Pee) and conversion probabilities (P^ 
evaluated using the Runge-Kutta numerical integration. 



P 1 



are 



^We use the standard parameterization [46] for the U matrix and the central value sinOi^ — 0.13 claimed 
in rcf. 1251. 



2.2 Neutrino electron scattering detection rates 

For the detection in SuperKamiokande and SNO through z/q, e~ — )■ z/^ e~ scattering, the NSI 
information comes in the probabihties and the cross section 
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where gL,R are the g^ji couphngs modified according to [llj'^l 
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with a = e, /i, r. 

For z/g both charged and neutral currents are possible, so that 
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whereas z/^ ,- only interact through neutral currents, hence 
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These expressions are then inserted in the spectral event rate 
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(17) 
which will be evaluated in the next section. Here (f){E) denotes the neutrino fiux from Boron 
and hep neutrinos, f{E^, E^) is the energy resolution function for SuperKamiokande and 
SNO [IZlllH] and the rest of the notation is standard. 

Notice that whereas the Hamiltonian (fT3!l is symmetric under the interchange 

^a/3 ^ ^f/3 for e, M, d 

such is not the case for the detection process [see eqs. (fT^ - (IT6|) ]. We finally note that at 
the detection level the NSI couplings e^^ are considered separately, as clearly seen from 
eq. (fT4l) . whereas at the level of propagation, since the diagrams involved in the interaction 
potentials add up, their sum should instead be considered. 



^This modification is different from the one in ref. [TT] since each £ coupling as defined in section [2] 
pertains not to a single vertex but to a whole NSI diagram. 



3 NSI couplings, probabilities and spectra 

We now perform an investigation of the effect of the NSI couphngs e^'o' = |£^'^' |e "Z' on 
the neutrino probabihty and event rates. Our aim in this section is to find those couphngs 
which lead to a flat SuperKamiokande spectral rate, thus improving the fit with respect 
to its LMA prediction while keeping the quality of the other solar event rate fits. We 
start by inserting one coupling at a time, first with equal CC and NC couplings, namely 
(£qS)cc = (£«pAfc = £«^ (subsection ETIl) and next with {e^J'p)cc ^ {£a^^)NC (subsection 
13.21) . We then extend our analysis to include all diagonal couplings in subsection 13.31 It is 
found that complex diagonal entries in the NSI Hamiltonian are the only ones the modify 
the LMA solution. 

For the sake of clarity we will organize the NSI couplings in three matrices according to 
whether the charged fermion in the external line is e, m, d 



(1^ 



Each set of three couplings e^'^' enters in equation flTTj) in the entry Vai3 of the interaction 
potential matrix. Altogether there are 18 couplings with 36 parameters: each matrix of 
the three in eq. flTSj) contains 6 independent entries, each with a modulus and a phase. We 
analyse one coupling at a time, by taking all others zero. We first consider the cases of 
purely real and imaginary couplings, hence 4 possibilities for each phase 




Le,u,d 
^a(3 



0, 7r/2, TT, (3/2)7r. (19) 



Motivated by the arguments expound in the introduction we investigate the parameter range 
|e„/3| e [5 X 10"^ , 5 x 10"^]. We find that 

• Off diagonal entries e^'^' (a 7^ (3) which contain 3x3x4=36 possibilities for moduli 
and phases do not induce any change in the LMA probability, nor any visible change 
in the rates, either if one or more at a time are inserted. 

• Diagonal entries e'^^''^ ^ . 

(a) Real couplings e'^^^''^ ^ = ilej^'^'*^ ^| (3x3x2=18 possibilities) do not change the 
LMA probability, hence the rates. 

(b) Imaginary couplings ^q^''^ ^ = ±i\e'^^''^ ^\ (3x3x2=18 possibilities) lead to prob- 
abilities which diverge from Plma for all \eaa\ > 5 x 10^^. According to the 
probability shape that is obtained, we group these cases in the following way 
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Table 1: The NSI couplings that modify the LMA probability. 



All cases in the first column of table [H along with D2, E2, F2 in the second column lead 
qualitatively to the same monotonically decreasing probability curve: a high probability 
(P > Plma) for low energy {E <Z MeV) and a low one (P < Plma) for intermediate and 
high energies. The curve becomes increasingly fiat in this energy sector as Saa increases, 
which is also reflected in the flatness of the SuperKamiokande spectral rate. However for 
cases Al, Bl, CI and D2, E2, F2 the probability gets too high for low energies so that the 
Ga [iniEO] rate fails to be conveniently fitted. The 'best' results in the sense that they lead 
to the most flat spectral rate which approaches the SuperKamiokande one and to a correct 
fit for Ga are obtained alternatively from cases Dl, El, or Fl for the following values 
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-i 1.5 X 10-3 
-i 2.5 X 10-3 
+i 2.0 X 10-3. 



(20) 



For larger values of the NSI couplings the probability moves further away from its LMA 
profile so that the Ga rate becomes too high and the ^B one too low. In figJUwe plot several 
survival probabilities: the dashed line is the vacuum one, then at the lowest energy and 
from bottom to top the first curve is the LMA one, the next corresponds to all cases in ( 120|) 
and leads to the best fit of the four, the next one to the case e^^ = —i 3x 10-^ and the top 
one to el/ = +i 5 x lO'^. A comparison is shown in table [2] between the predictions and 
the quality of the fits obtained from LMA and the case e'l^ = —i 1.5 x lO"^. In figsj2]and 
|3] we show the SuperKamiokande spectrum for LMA (upper curves) and for the first case 
in fl20|) (lower curves) superimposed on the data points taken respectively from refs. [29] 
and [30]. The improvement obtained through the NSI coupling is clearly visible. In figJH 
the two curves are superimposed on the SNO data points for electron scattering [5l]. Here 
the data are also clearly consistent with a constant rate. 

The first set of cases in the second column of table [Tj namely A2, B2, C2, lead qualita- 
tively to the inverse behaviour with energy of the LMA probability. As \eaa\ increases from 
its lower bound, one gets P < Plma for low energies (P < 2 — 3 MeV) and P > Plma for 
intermediate and high energies, so that the fits worsen with respect to the LMA ones. 

Furthermore all cases in the third column of table [H namely A3, B3, C3, D3, E3, F3 
lead to probability curves which are totally unsuitable: they deviate drastically from both 
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Table 2: Comparison between the LMA predictions for solar event rates and the NSI ones 
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with — ^l^^^/'l = —i 1.5 X 10 ^. For details of the x^ analysis see for instance [211/ . 



Plma and a from fiat, suitable profile able to generate the SuperKamiokande spectrum. 

We have also checked that combinations of real and imaginary parts for all couplings do 
not change the previous results. This should be expected since, as mentioned earlier, purely 
real couplings do not change the LMA probability. The only consequence of introducing 
real parts in the NSI couplings comes in the spectral event rates through the quantities gL,R 
in eqs. flT^ and flT5]) . but the differences lie much beyond the experimental accuracy. 

The analysis of the more general case of different CC and NC couplings affecting the z/g e 
scattering diagrams can be made quite simple if one examines the coefficients of {s^ap)cc 
and {£^ap)NC in eq.(ll). The first is unity whereas for equal CC and NC couplings it is 0.96 
and the second is now -0.04 as compared to the previous value 0.96 as well. Consequently 
one expects that the analysis for {£^ap)cc leads to approximately the same results as for 
equal couplings while the results are modified by a factor of 0.96/(-0.04) in the analysis for 
(£^^)7vc. Indeed the convenient modification in the LMA probability is obtained for 



or alternatively 



K^)cc = -^{K,''\)cc = -^ 1-4 x 10^^ (21) 

K^)nc = +^(l</l)7vc = +^ 3.6 X 10-2. (22) 

which lead to the same probability as the cases listed in eq. fl20|) . The results for the other 
couplings involving u and d quarks are of course unchanged. As before we have considered 
one coupling at a time to be non zero. 

3.3 All diagonal couplings 

We now allow for all diagonal entries of the NSI term in eq. (IT2|) to be nonzero and complex. 
In a first step to this generalization we start with all three couplings e!^'^''^ nonzero simulta- 
neously. Starting with the three equal in moduli, we find that the right modification of the 
LMA probability is obtained with 

^^/ = -^k^/l = -^0.7x10-3 

^m/ = -^k;:/l = -^0.7x10-3 (23) 

4/ = +^I4/I = +^ 0.7 X 10-3 _ 



In a second step we allow for one diagonal entry in the NSI term of (fT2|) in addition 
to sf^J^''^ to be nonzero and complex. We find that a finite e^'"'*^ as added to the parameter 
choice ( !23l) leads to a 'wrong' probability, unless |e^'"''^| < |(1/100)£:^'^''^| in which case the 
LMA solution remains unchanged. The same is true for el'^'"^ 7^ 0. However adding both 
^e^M,d ^^^ e^'^''^ to (125]) . the correct change in the LMA probability can be obtained provided 
\el'^''^\ and le^'^'"^! ^ |(l/10)e!^'"''^|. To this end two choices are possible: either the signs of 
g^,u,d g^j,g changed with respect to those of table [1] with those of e^'"''^ unchanged that is 

f"" = -yV I r" = -?V" I f"^ = -\-7\f^ I ('24') 

and 

4, = -z|e^J, e:^, = -^ICI, et = +^k'.l (25) 

or vice-versa. 

A solution with all nine diagonal couplings having equal moduli is also possible. It 
corresponds to 

kef '1 - I^m'I - Kr '1 = (2 - 4) X 10^^ (26) 

provided the signs of fI'^''^ are unchanged with respect to table [1] and those of e^'"'"^ are 
reversed that is 

f-e ^ +?V I f" = +?V I f'^ = — ?V^ I ('27') 

and 

Notice that the reverse choice in signs with respect to fl7r|) and f l25]) would leave the LMA 
solution unchanged. 

In this section we have analysed the conditions to be imposed on the NSI couplings 
in order to obtain a suitable survival probability for a flat SuperKamiokande spectrum, an 
improved CI rate prediction relative to the LMA one and accurate predictions for all other 
rates. As shall be seen next, these conditions while necessary, are not sufficient, owing to the 
fact that neutrino decay in solar matter follows as a consequence of the imaginary diagonal 
couplings in the Hamiltonian. 

4 Neutrino decay in solar matter? 

We recall that a stable stationary state solution of the wave function, for a particle of 
energy E, contains a phase factor e~*^*. If the Hamiltonian has complex eigenvalues this 
phase factor becomes instead 

so the state is unstable, with decay rate F and lifetime F"^. 

Since the suitable Hamiltonian that reduces the tension between the LMA solution and 
the data has imaginary diagonal couplings, its eigenvalues are complex which implies the 



existence of unstable states with a finite decay rate. We note that this instabihty is induced 
by the solar matter density [see eq. (TT2|) ] and in its absence, in the vacuum, the neutrinos 
are stable. Furthermore, as the number of neutrinos (or neutrinos and antineutrinos) must 
remain constant as a consequence of unitarity, we must look for those solutions of the wave 
equation in which at least one but no more than two eigenvalues have a negative imaginary 
part. The remainder will have positive imaginary parts. This will ensure that the unstable 
state or states decay into any of the lighter. The negative imaginary parts correspond to the 
decaying states as it follows from f l29l) . while the positive ones to the states into which the 
former decay. Such criteria must be imposed together with the requirement of consistency 
with the solar neutrino data, namely a fiat SK spectrum and accurate Ga, CI and SNO 
rates. To this end the correct solution is found to be 



Hmsx = GfV2N, 



-F \ I -F \ I — -F 

2 \ / 2 \ / 2 



where f = 3.5 x 10 . Here 



xi I -tF +^2 ~^£ "^^3 ^^ 
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(30) 



1 2n 8 2n Ar„ I 2 ^^ Nn , ^ 

a^i = 2 + 2sm Ow, ^2 = 1- -s%n 9w + ^j^, X3 = -- + -sm ^w - ^ (31) 
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Notice that ([30]) is consistent with (|26D. fl271l and fl28|l . 

In figj5] (panels (a) and (b)) we plot the real and imaginary parts of the eigenvalues 
of the Hamiltonian flT2|) with the solution fl30|) for the NSI term as a function of the solar 
fractional radius for a typical solar neutrino energy E = 1 MeV. From panel (b) which 
displays the imaginary parts, it is seen that a significant decay of one of the mass matter 
eigenstates can occur during the first 30-40% of the trajectory with an average decay rate 
O(10~^^ eV). Further along the neutrinos appear to be stable while still in the sun and in 
the vacuum. As expected, the decaying state (the negative values in panel (b)) corresponds 
to the largest of the mass matter eigenstates (panel (a)). 

The neutrino decay modes considered in the literature involve photon or majoron emis- 
sion along with a neutrino or antineutrino. Although radiative decay is enhanced in a 
medium [52], for the sun such an enhancement is far too small, as it leads to a lifetime 
at least five orders of magnitude larger than the age of the universe, so such an effect is 
totally irrelevant here. The same holds for the 'neutrino spin-light' in matter [53l|55. We 
are therefore left with the possibility of neutrino decay in matter into an antineutrino and 
a majoron x [551458] 

Ui ^Uj + X (32) 

where subscripts i and j denote mass eigenstates. 

In the energy range of interest within the sun the neutrino transition is not matter 
dominated, hence mass and flavour eigenstates do not coincide. So, having started from the 
interaction potentials which are flavoured based, a rotation from flavour to mass states in 
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solar matter is needed in order to obtain the decay rate. However the quantity of interest 
which can be directly related with experiment is the antineutrino probability at the Earth 
for a particular flavour, namely Ve- To this end the above mentioned rotation must be 
undone in the vacuum. Thus we will work from now on with flavour based quantities. 

The probability of antineutrino appearance for flavour /3 per unit solar radius is then 
'^i^^ = f"°- ME«) (1 - e-('.-.)') f!£(l#j£) ,E„ . (33) 



dr JEf dEf 



Here 0(-Eo) is the normalized ^B solar neutrino spectrum [59] with Eq^^^ = 16.56 MeV and 
the factor (1 — e"'"'^^'-^-'')'') takes into account the neutrino flux reduction. The quantity -^ 
denotes the differential decay rate per unit antineutrino energy Ef 



dT _ \gai3\ Eq — E 
dE'f ~ "Svr E^ 



f 



2 



V^ir)-Vf,ir)\{NSI) (34) 



where gai3 are the neutrino- majoron couplings and V^, VjS the neutrino and antineutrino 
interaction potentialqj- The energy dependent antineutrino production probability from the 
decay ( 132|) is therefore 

^..(^/) = ^/ \VAr)-y,{r)\{NSI)\ <^{E,){\-e-^^^^^-y)^-^ dE.dr (35) 

^T^ jRi JEf -C/Q 

where Rs is the solar radius and Ri is the production point, assumed for simplicity at 5% 
of Rs for ^B neutrinos. Upper bounds for the neutrino majoron couplings gaj3 have been 
estimated [55l[58] and in particular the authors of ref. [58] flnd 

^l^eaP <5.5x 10"^ (36) 

a 

of particular interest here, since we will be concerned with z/g production for comparison 
with the Borexino [HD] and KamLAND [HI] upper bounds on solar v^- The probability 
(|35|) is calculated in the interval [57] [max{l/2(V^ — V p)^ ~^ p}^ ^OmaA ^^^ the decay f l32|) 
requires ImiVa — V p) > with Va, V p obtained from eqs. fITT]) . f l30|) and fl3T]) 

V,{NSI) = '-GpV2{3.5xlO-^)Njl-^stn'ew + ^) (37) 

V^iNSI) = tGFV2{3.5xlO-'')Nj-l + Pin^ew-^] (38) 

KiNSI) = '-GpV2{3.5 X 10-^)iV, (l - ^-sm^Q^ + ^"j . (39) 



^Recall that neutrino and antineutrino interaction potentials are symmetric to each other, V/3 = —Vj^ 
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We next evaluate the interaction potentials (l37I) - (l39l) throughout the neutrino trajectory in 
the sun (see figj6]), and analyse the possibilities for neutrino-antineutrino decay in terms of 
the necessary condition Im{Va — V13) > 0. To this end, inspection of figj6] shows that the 
sums Ve — Ve, Vr — Ve and V^ — Ve (the two top curves and the fourth from the top in panel 
(a)) contribute along the whole trajectory to Ve production. Similarly v^ production is less 
important, as it can depend only on Vg — V ^ (fourth line from the top in (a)). Finally z/^ 
production comes from the contribution of K — Vr and V-r — Vr (second and third curves 
from the top in panel (a)) which are both positive over the whole trajectory. The two bottom 
curves in panel (a) namely V^ — Vr (or 1/^ — Vfj) and V^ — V^ are both negative, they refer 
to Vr and z/^ production and need not concern us here. To summarize, the production of: 

• Ve involves Ve — Ve, Vr — Ve and V^ — Ve 

• v^ involves Ve — V^ 

• Vr involves Ve — Vr and Vr — Vr- 

Following these rules and using eq.( l35|) with the upper limit (136|) . we represent in figE] 
the upper bound of the energy dependent Ve production probability, Pp^{Ef). In panel (a) 
we use a logarithmic scale and in panel (b) both the LMA probability and Pp^{Ef) (inner 
graph) are shown. The antineutrino probability is seen to be extremely small but grows 
rapidly as Ef approaches its lower limit. It decrease fast to zero as Ef — )• Eq,^^^, since fewer 
neutrinos contribute in the upper energy range. 

In order to seek a comparison with the data, we first note that the Borexino upper 
bound on solar Ve, namely 760 cm~'^s~^, applies for Ep^^ > 1.8 MeV which includes most of 
the ^B spectrum ^60j and so corresponds in practice to a maximum total probability for z/g 
production [ID] 

^- ^ i3^ - '■' - ''-'■ («) 

We therefore consider the area limited by Pp^{Ef) and the two extreme abscissa 1.8 MeV 
and 16.56 MeV, dividing it by the area limited by the unit probability and the same two 
abscissas. This ratio gives the total probabihty for Ve production to be compared with fj40|) . 
We find for the ratio of these two areas 

7.5 X 10"^^ (41) 

which is within the Borexino bound by 6 orders of magnitude.. 

The KamLAND collaboration, on the other hand, reports a bound 370 cm~^s~^ for 
Ep^ > 8.3 MeV [BT] which corresponds to a probability 

370 

Pu < Tw^ = 2.1 X 10"^ (42) 

^"5.94x106 J^'^'' <PiE)dE ^ ' 
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whereas for the ratio of the corresponding areas we now get 

3.3 X 10^12 (43) 

within the KamLAND bound by 7-8 orders of magnitude. Hence any increase in experi- 
mental sensitivity will be unable to reveal a possible solar antineutrino flux produced from 

NSI. 

Finally the full physical process in our model for neutrino propagation and decay through 
NSI in the sun is represented in figJHl Equation (TT^ we used for the Hamiltonian does not 
take into account the extra physics involved in the majoron coupling and is therefore a 
truncated Hamiltonian, whose hermiticity is restored once the detailed majoron emission 
process is taken into account. 



5 Conclusions 

We have investigated the prospects for improving the LMA predictions for solar neutrino 
event rates with NSI. At present there is no evidence of any new physics associated to a scale 
not too far above the electroweak scale, hence the great variety of theoretical models available 
for NSI. In our approach we assumed that NSI are extra contributions to the vertices z/q,z//3 
and z/qC, so the new couplings describe the deviation from the standard model. With this in 
mind we derived the neutrino interaction potential in solar matter which was added to the 
standard Hamiltonian, proceeding with the integration of the evolution equation through the 
Runge-Kutta method. Neutral and charged current couplings are involved in interactions 
with electrons whereas only neutral couplings affect those with quarks. We considered the 
new interactions both at the propagation and at the detection level. The improvement we 
searched for the LMA predictions consisted in finding whether and how the modification 
induced by NSI can lead to a fiat spectral event rate for SuperKamiokande and an event 
rate for the CI experiment within la of the data, while keeping the accuracy of all other 
predictions. 

We used the current notation for the NSI couplings e^'^' where a, /3 are the neutrino 
labels, e, u, d denote the charged fermion involved in the process and we investigated the 
range |£:q,^| e [5 x 10~^ , 5 x 10~^]. Our most remarkable and intriguing finding is that 
only the imaginary diagonal couplings in the Hamiltonian provide a solution to the tension 
between the LMA prediction for the SuperKamiokande spectrum and the data. This implies 
the existence of unstable neutrino states, the instability being induced by the solar matter. 
In the vacuum the neutrinos remain stable. 

We may now summarize our main results as follows 

• Diagonal, imaginary couplings e^'o''^ ^ = ±i\e^'^''^ ^\ are the only ones that lead to 
changes of all kinds in the LMA probability and hence the rates. 
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Real couplings e 



e,u,d P 

al3 



±\e 



e,u,d P\ 
a/3 



do not change the LMA probability when consid- 



ered either one at a time or altogether and thus they induce a small change in the 
neutrino electron scattering rate (< 1%) which is far beyond experimental visibility. 
For the same reason the addition of real parts to the imaginary couplings does not 
change the results. 

Off diagonal couplings e^'^' {a ^ (3) whether real or imaginary considered either 
one at a time or altogether do not change the LMA probability and thus the rates in 
a significant way. 

Solving the tension between the LMA solution and the data, in particular predict- 
ing a flat SuperKamiokande spectrum, requires unstable neutrinos, owing to the fact 
that only imaginary diagonal couplings in the NSI Hamiltonian modify LMA. The 
requirement of unitarity then gives the following solution 



e,u . 
ee 
e,u . 

MM ' 

e,u . 

TT 



-IE 



"ee 

-d 

''ii.il. 



2"- 
%E 



(44) 



with e = 3.5 X 10" 



• The above facts imply neutrino decay from the heavier to the lighter states in the 
solar matter probably with majoron emission, since radiative decay is irrelevant in the 
sun. Having calculated the antineutrino appearance probability, we flnd however that 
it is quite unlikely to ever observe experimentally antineutrinos from the sun due to 
NSI. Small diagonal imaginary NSI couplings lead to quite a small and unobservable 
antineutrino production probability. However this induces a visible and remarkable 
change in both the survival and conversion probabilities to z/^ and z/,-. 

• Experimentally the flatness of the SuperKamiokande and SNO electron energy spectra 
is at present the only evidence in favour of neutrino NSI in the sun. 

Finally the detailed physics involved in majoron emission was not taken into account 
and the Hamiltonian we used eq.( ll2p is therefore a truncated one which is sufficient for our 
purpose. 
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Figure 1: Survival probabilities as a function of neutrino energy in eV. The dashed line is 
the vacuum one. The remainder are (a) the LMA one, the one providing the best fit to the 
data (b) with e^^^ = —i 1.5 x 10~^ [eg. ^2^) 1. curve (c) fore^^^ = —i 3 x 10^'^ and curve (d) 
for e^/" = +i 5 X 10"^ with other non standard couplings vanishing in each case. The last 
two curves ((c) and (d)) lead to an unacceptably high Ga rate prediction ^5l[5Z 
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Figure 2: Predictions for SuperKamiokande (units in MeV for electron energy). The upper 
curve is the LMA spectrum and the lower curve is the LMA spectrum with non standard 
interactions as in eqs.^2^ or (24^ . These are superimposed on the data published by the 
Collaboration in 2002 !M 



20 



B 

DC 

■D 
0) 
O 

■D 
CD 

DC 

::^ 

CO 



0.62 r 

0.58 'r 




0.34 t 
0.30 



10 15 

Ee (MeV) 



20 



Figure 3: Same as figlE with the data published in 2008 |S 
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Figure 4: Predictions for SNO neutrino electron scattering superimposed on the data p7]/ 
(units in MeV for electron kinetic energy). The upper curve is the LMA prediction and the 
lower curve is the LMA one with non standard interactions. Error bars are larger than in 
SuperKamiokande so that the data are consistent with a flat spectrum. 
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Figure 5: The real (a) and imaginary parts (b) of the neutrino mass matter eigenvalues for 
E = 1 MeV : the lower curve in (b) associated to the decaying state corresponds to the upper 
curve in (a), the largest of the mass matter eigenvalues. 
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Figure 6: Neutrino NSI potentials along the sun: panel (a) from top to bottom displays the 
imaginaryparts of K - V^, K - V^ (or K - Ve), Vr-V^, K - V^ (or V^ - V^), V^ - V^ 
(or Vr — Vfj,) and V^ — V^^. The first, second and fourth from the top, pertaining to v^ 
production, are also shown in panel (b) in logarithmic scale. 
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Figure 7: Panel (a): antineutrino production probability in a logarithmic scale. Panel (b): 
LM A +N SI probability and antineutrino production probability (inner graph). Its extremely 
small value prevents the possibility of observing antineutrinos from the sun from NSI. 

25 



V. 



V 



e,iLij 



V 



e,iij 



W 







V 



e,/i,r 



CyUyd 




e,u,d 



Figure 8: The processes involved in the propagation of neutrinos in the sun as in our model: 
the two upper diagrams are the standard ones for matter oscillation and the lower one 
represents the decay z/j — )• Vj + majoron{x)- 
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